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The molecular dynamics study of fragmentation in peripheral ^®^Au +^^^Au collisions at 35 
MeV/nucleon is presented to probe the nuclear matter compressibility in low density regime. The 
yields of different fragment species, rapidity spectra, and multiplicities of charged particles with 
charge 3 < < 80 are analyzed at different peripheral geometries employing a soft and a hard 
I equations of state. Fragment productions is found to be quite insensitive towards the choice of 

, nucleon-nucleon cross sections allowing us to constrain nuclear matter compressibility. Comparison 

_ ' of calculated charged particle multiplicities with the experimental data indicates preference for the 

soft nature of nuclear matter. 
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^ ' I. INTRODUCTION 

During the last two decades, numerical simulations of medium and high energy heavy-ion (HI) reactions have 
provided a unique opportunity to explore the nuclear matter at the extreme conditions of density and temperature 
! UHll ■ The compression of nuclear matter can be judged via equation of state (EoS) which is also a main input for 
any theoretical model along with nucleon-nucleon {n-n) cross section 04^- Various attempts have been made to find 
the observables which are sensitive to the nuclear EoS. In the past, comparison of theoretical predictions with the 
. experimental results has been used to extract the nuclear equation of state. One of the earlier attempts for nuclear 
' EoS (or, incompressibility) was via giant monopole resonance (GMR) studies The scattering of a particles off 

: the nucleus induces volume oscillations with L=0, which can be used to determine the incompressibility 'k' of that 
(<~^ ' nucleus. These results generally yield incompressibility in the range k ~ 250-270 MeV indicating the matter to be 
softer. A recent GMR study in the ^°^Pb and ^°Zn nuclei showed that the softening of nuclear matter is needed to 
J* • explain the collective modes with different neutron-to-proton ratios . Another study on the fusion reported linear 
momentum transfer to be sensitive to both the EoS and n-n cross section (l2j . Within quantum molecular dynamics 
(QMD) model, an incompressibility of k=200 MeV (i.e.. soft EoS) was reported to reproduce the experimental data 
■ on energy transfer in a compound nucleus formation |12| . 



Collective fiow observed in HI collisions is another observable which is found to be sensitive towards the stiffness of 
nuclear EoS J^, H, [l^ . The collective transverse in-plane fiow and balance energy (the energy at which fiow becomes 
zero) have been studied extensively over the past two decades so as to constrain the EoS, but still the uncertainties 
are very large. For example, a stiff EoS with k=380 MeV regroduces the transverse flow data equally well as obtained 
with soft momentum dependent EoS with k—210 MeV jsl. Similarly, comparison of transport model calculations 
with data of EOS Collaboration for the energy dependence of collective flow favored neither the 'soft' nor the 'hard' 
equation of state Q- In recent comparison of elliptic flow data with microscopic transport model calculations of the 
Refs. no consistent agreement with the data could be obtained [17] for the two different models of Refs. [15| 



& [16[ . Therefore, it is clear from the above review that an appropriate choice of the nuclear equation of state is still 
far from settlement. The task of deriving quantitative information about the EoS requires detailed comparison of 
theoretical calculations assuming different equations of state with experimental data. 
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At lower beam energies, Pauli blocking of final states gets more pronounced. As a result, mean-field effects and long 
range Coulomb force govern the reaction dynamics. If one goes to still lower energy regime, well-known phenomena 
like complete fusion, incomplete fusion, fission, cluster emission etc., can be seen At incident energies above 

20 AMeV, phenomena like production of intermediate mass fragments (IMFs), projectile-like and target-like fragments 
(PLFs and TLFs) dominate the exit channel. The phenomenon of multifragment-emission in low energy domain is, 
however, least exploited to infer the nuclear EoS. Naturally the study of fragment-emission in low energy domain may 
be of importance to constrain the nuclear incompressibility, where the role of different n-n cross sections is expected 
to be minimal. To explore the possibility of achieving information on the nuclear EoS, we plan here to simulate 
the peripheral reactions of ^^''Au -|-^^^Au at Eiab^35 AMeV and at different peripheral geometries where accurate 
data has been measured on Multics-Miniball set-up j20|. To this end, we performed detailed calculations within the 
quantum molecular dynamics (QMD) model [2ll [22j which is described in detail in section |TT] along with simulated 
annealing clusterization algorithm (23l . [23 | . Section IIIII presents results of our numerical calculations and comparison 
with available experimental data, which are finally summarized in section ITVl 



II. THE MODEL 



A. Quantum Molecular Dynamics (QMD) Model 

The quantum molecular dynamics model is a n-body transport theory that incorporates the quantum features of 
Pauli blocking and stochastic n-n scattering. Each nucleon in the colliding system is represented by a Gaussian wave 
packet as: 



^,(r,p,(t),r,(i)) = ^27,-L)3/4 



(1) 



Mean position ri{t) and mean momentum Pi{t) are the two time dependent parameters. The Gaussian width \/L is 
fixed with a value of 1.8 fm and is same for all nucleons. This value of \/L corresponds to a root-mean-square radius 
of each nucleon. The centers of these Gaussian wave packets in TZ^ and Vs spaces follow the trajectories according to 
the classical equations of motion: 

The Hamiltonian T-l appearing in Eq. ([2]) has contribution from the local Skyrme-type, Yukawa and effective Coulomb 
interactions [2ll |: 

V^°' = tiSir,-r,) + t2S{r,-r,)6{r,-rk) (3a) 
, ^ ^^eM-]r^-r^ ^^^^ 



Zi ■ Z i e 
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= -1 — 1^. (3c) 
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Here Zi, Zj are the effective charge of baryons i and j. The long-range Yukawa force is necessary to improve the 
surface properties of the interaction. The parameters /z,ii,i2,^3 appearing in Eqs. ([3]) are given in Ref. [2^. These 
parameters are adjusted and fitted so as to achieve the correct binding energy and root mean square values of the 
radius of the nucleus [l^. In QMD model, one neglects the isospin dependence of nucleon- nucleon interaction. All 
nucleons in a nucleus are assigned the effective charge Z = ^^^^^ [2l| . It is worth mentioning that isospin dependent 
flavor of QMD model {i.e. IQMD) has also been used in literature [ll], [2^. This microscopic transport code explicitly 
takes into account the differences of neutron and proton potentials and cross sections. The Skyrme part of the 
interaction used in QMD model has the generalized form: 

'PiA ^ P f Pij 
2 V Po / 7 + 1 V Po 



= ■ (4) 



The interaction density in Eq. (|4]) is defined as: 



(47rL)3/4 



g-(r.-r,)V4L^ (5) 
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The parameters a and /3 are adjusted to reproduce the infinite nuclear matter binding energy (E/A = -16 MeV) at 
saturation nuclear matter density p^. The third parameter 7 can be varied independently to account for different 
nuclear incompressibilities k (i.e. different equations of state). Two different parameterizations are used: a soft EoS 
with incompressibility k=200 MeV, and a hard EoS with k=380 MeV. The parameters a, /3, and 7 employed in the 
QMD model serve as an input for repulsive potential with high compressibility (i.e. hard EoS), and less repulsive 
potential (i.e. soft EoS). The standard parameters corresponding to these two equations of state are listed in TablelH 



TABLE I: The Skyrme parameters for 'soft' and 'hard' interactions used in the QMD model. 



EoS Q(MeV) /3(MeV) 7 ^(MeV) 
Soft -356 303 7/6 200 
Hard -124 70.5 2 380 



The influence of different n-n scattering cross sections will be determined by employing a set of different cross 
sections varying from energy-dependent cross section [26| to constant and isotropic cross sections with magnitudes of 
40 and 55 mb. It is worth mentioning that in recent times, even relativistic version has also been analyzed [l!|. As 
noted in Ref. 1], this has no effect on present findings. A hard EoS with energy dependent cross section is labeled as 
HardP^ . Incorporation of isotropic and constant cross sections of 40 and 55 mb strengths have been labeled as Hard'^'^ 
and Hard^^ , respectively. Similarly, for the soft equation of state, we have Soft^^, Soff^^ and Soft^^, respectively. 
Since QMD model follows the time evolution of nucleons only, one has to employ secondary clusterization algorithms 
to identify fragments' structure. In the present paper, simulated annealing clusterization algorithm (SACA) has 
been used to identify final fragment structure. This method is reported to explain the ALADiN data on spectator 
fragmentation quite nicely at relativistic bombarding energies [23 |. 



B. The SACA Formalism 



This clusterization procedure allows early identification of fragments before these are well separated in coordinate 
space. In the SACA method, fragments are constructed based on the energy correlations. It works on the principle of 
energy minimization of fragmenting system. Standard minimum spanning tree (MST) procedure [2^ is employed to 
obtain pre-clusters. This procedure assumes that nucleon pairs separated by distance | — |< 4 fm, belong to the 
same fragment [sl, [2^. Such preliminary cluster configuration is determined at every time step. Thus one can address 
the time evolution of mass, charge, position and momentum of each fragment or single nucleon. The pre-clusters 
obtained with the MST method [22, [l^ are subjected to a binding energy check [23| : 



^ rv — 1 



Nf 
1 \ 

^ Van (ra-r^) 



hind-) 



(6) 



with E] 



hind 



4.0 MeV if iV/ > 3 and 



^bmd = otherwise. In Eq. (|6]), Nf is the number of nucleons in a fragment 
and PJv™' is the center-of-mass momentum of the fragment. The requirement of a minimum binding energy excludes 
the loosely bound fragments which will decay at later stage. To look for the most bound configuration (MBC), we 
start from a random configuration which is chosen by dividing whole system into few fragments. The energy of each 
cluster is calculated by summing over all the nucleons present in that cluster using Eq. ([6]). 

Let the total energy of a configuration k be Ek{= J^i^fd)' where Nf is the number of nucleons in a fragment and 
is the energy per nucleon of that fragment. Suppose a new configuration k (which is obtained by (a)transferring 
a nucleon from randomly chosen fragment to another fragment or by (b) setting a nucleon free, or by (c) absorbing 
a free nucleon into a fragment) has a total energy E^,/ . If the difference between the old and new configuration 
A£^(= E/.' —Ek) is negative, the new configuration is always accepted. If not, the new configuration k may nevertheless 
be accepted with a probability of exp{—AE/v), where v is called the control parameter. This procedure is known 
as Metropolis algorithm. The control parameter is decreased in small steps. This algorithm will yield eventually the 
most bound configuration (MBC). Since this combination of a Metropolis algorithm with slowly decreasing control 
parameter v is known as simulated annealing., so our approach is dubbed as simulated annealing clusterization algorithm 
(SACA). We have used here an extended version of SACA, in which each cluster is subjected to its true binding energy 
based upon modified Bethe-Weizsacker mass formula [28[. It may be stated that fragmentation analysis performed 
within extended version yields the same results as with constant binding energy check of -4 MeV/nucleon. It also 
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justifies using extended approach to analyze the fragmentation at such a low incident energy when nuclei are still 
in Fermi energy domain. The constant binding energy criterion (of -4 MeV/nucleon) is chosen keeping in mind the 
average binding energy of clusters. We have also analyzed the fragmentation pattern employing Ebind based upon 
experimental binding energies. Nearly no effect of this modification was found. For further details, we refer the reader 
toRefs. [2l[23,|23,[2|. 

III. RESULTS AND DISCUSSION 

Figure 1 shows the time evolution of ^^^Au +^^^Au reaction at 35 AMeV with reduced impact parameters 
b/bmax=0.55 (left panel) and 0.85 (right panel). Top panel depicts the time evolution of average nucleon density 
pavg fQj. gQf^ g^j^^ hard EoS. One notices several interesting results: (i) Maximal density is reached nearly at the same 
time at both impact parameters, whereas saturated values is slightly more for higher impact parameters, (ii) Choice 
of different n-n cross sections have insignificant influence on the results obtained. This happens due to effective Pauli 
blocking at such a low incident energy that prohibits n-n collisions. 

Stiffness of nuclear EoS, however significantly influences the mean nucleon density and other fragment observ- 
ables shown in subsequent windows. This difference is clearly visible in the evolution of heaviest fragment (yl™"^), 
multiplicities of free particles, light charged particles LCPs [2 < A < 4], and clusters with mass A> 5. The mean size 
of heaviest fragment (A™°^) attains minimum around 100 fm/c, where stable fragment configuration can be realized 
and compared with experimental results. With stiff EoS, heavier (^™'^^) is registered. Interestingly, multiplicity of 
free particles obtained also follows the same trend as (A™"^^). This means that dissipation of energy takes place 
mainly via emission of free-nucleons that cools down the nuclear system in case of hard EoS. Consequently, lesser 
yields of LCPs and fragments with mass A> 5 are obtained with a stiff EoS. On other hand, soft EoS favors emission 
of LCPs and heavier fragments {A > 5) from spectator zone, thereby decreasing the size of A"^°-^ . The insensitivity of 
fragmentation pattern towards choice of different n-n cross sections may be, therefore, useful to constrain the nuclear 
compressibility in low density regime. 

Next we study the rapidity spectra of free nucleons and intermediate mass fragments in transverse and longitudinal 
direction using a hard and a soft equations of state. The transfer of excitation energy from the participant zone to 
spectator matter has a direct bearing on the rapidity distribution of the fragments. Figure 2 displays the spectrum of 
scaled transverse y^^-* and longitudinal y'^^' rapidity distribution of free particles (top) and intermediate mass fragments 
IMFs [h < A < 65] (bottom) for the collision of Au (35 AMeV)+ Au at reduced impact parameter h/hmax=^-^^- 
As expected, the rapidity spectrum of free nucleons and IMFs is quite sensitive to the nuclear EoS that brings out 
significant change in their transverse expansion as well as stopping pattern. Using hard interactions, a larger fraction 
of free nucleons are emitted into transverse direction. IMFs aren't, however, dispersed much into transverse directions 
and continue to move at target and projectile velocities. Similar trends are visible for longitudinal rapidity (j/z) 
distribution as well. Using a 'stiff' EoS, system seems to cool-off via abundant production of the free nucleons from 
the midrapidity as well as from the spectator zone, whereas a 'soft' EoS contributes significantly towards the IMFs 
emission at target and projectile rapidities. It means that system propagating under the soft interactions is less 
equilibrated. As a result, heavier fragments leave the participant zone quite early and suffer less collisions. These 
findings suggest that fragment emission from the decay of spectator component is quite sensitive to the mean field 
and compressibility of participant matter. 

Finally we calculate the multiplicity of charged particles with 3 < Z < 80 using a hard and a soft equations of 
state for Au(35 AMeV)-|-Au collisions at six peripheral geometries (see Fig. 3). The multiplicities calculated at 
100 fm/c are subjected to forward rapidity condition (y > 0.5 ybeam) m the center-of-mass frame to exclude events 
from midrapidity and quasitarget decay. Also shown in the figure is integrated multiplicities of charged particles 
with charge 3 < Z < 80 (i.e. J^^ N(Z)dZ) obtained on Multics-Miniball setup [20|. It is worth mentioning that 
multiplicities were calculated keeping in mind the angular range covered by combined Multics-Miniball array. Overall 
we can see that results obtained with a soft EoS are consistent with experimental data at all colliding geometries. 
These findings reflects the ability of molecular dynamics approaches (QMD, in our case) to describe the reaction 
dynamics in low-energy regime. 

Due to more explosive nature of hard EoS, spectator matter mainly de-excites via emission of free nucleons and 
therefore, decline in multiplicity of heavier clusters occurs. An increasing trend of fragment multiplicity with centrality 
can be understood in terms of more excitation energy deposited in spectator matter. In semi-peripheral events, a larger 
chunk of excitation energy gets transferred to spectator matter, thereby, leading to rise in multiplicity of fragments 
with decrease in impact parameter. Nuclear mean-field, therefore, becomes important factor governing the outcome 
of spectator decay, while nucleon-nucleon collisions dominate the participant matter physics. This analysis clearly 
illustrates the relatively soft nature of nuclear matter. 
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IV. SUMMARY 

In summary, QMD model has been used to infer the inter-play of different model inputs on fragment-emission in 

peripheral ^'^^Au +'^'''^Au collisions at 35 AMcV. Wc find that choice of different nuclcon-nuclcon cross sections has 
marginal role to play at such a low incident energy. However, the multiplicity of charged particles obtained from the 
spectator decay are strongly influenced by the incompressibility of the nuclear matter. The hard equation of state 
results in enhanced emission of free nucleons and fewer heavier fragments. Model calculations with soft EoS are 
found to give encouraging results which are in accord with experimental trends. This study favors soft nature of the 
nuclear matter. 
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Figure Captions 

FIG. 1. QMD simulation of Au (35 AMeV)+Au collisions at reduced impact parameter b/bmax—O-^^ (left panel) 
and b/bmax—0-&5 (right panel) as a function of time: (a) mean nucleon density p""^^ / Po\ (b) size of heaviest fragment 
j^max.^ multiplicities of (c) free nucleons, (d) light charged particles LCPs, and (e) fragments with mass A > b, 
respectively. 

FIG. 2. Rapidity distribution dN/dy of free nucleons and intermediate mass fragments (IMFS) as a function of 
scaled transverse, y'^^^ /yteam (left) and longitudinal, y^^^ /ybeam (right) rapidities in Au (35 AMeV)+Au reaction at 
reduced impact parameter 5/5,„aa;=0.55. Solid and dashed curves correspond to model calculations using a 'soft' and 
a 'hard' EoS respectively. 

FIG. 3. The impact parameter dependence of multiplicity of fragments with charge 3 < Z < 80 obtained using 
a 'soft' EoS (sohd line) and a 'hard' EoS (dashed line) in Au(35 AMeV)+Au collisions. Filled circles depict the 
experimental data points [20l |. 
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